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ABSTRACT NOMENCLATURE 


The wake behind a two-bladed model rotor in 
light climb was measured using particle image ve- 
locimetry, with particular emphasis on the develop- 
ment of the trailing vortex during the first revolution 
of the rotor. The distribution of vorticity was dis- 
tinguished from the slightly elliptical swirl pattern. 
Peculiar dynamics within the “void” region may ex- 
plain why the peak vorticity appeared to shift away 
from the center as the vortex aged, suggesting the 
onset of instability. The swirl and axial velocities 
(which reached 44 and 12 percent of the rotor-tip 
speed, respectively) were found to be asymmetric rel- 
ative to the vortex center. In particular, the axial 
flow was composed of two concentrated zones mov- 
ing in opposite directions. The radial distribution 
of the circulation rapidly increased in magnitude un- 
til reaching a point just beyond the core radius, after 
which the rate of growth decreased significantly. The 
core-radius circulation increased slightly with wake 
age, but the large-radius circulation appeared to re- 
main relatively constant. The radial distributions 
of swirl velocity and vorticity exhibit self-similar be- 
haviors, especially within the core. The diameter of 
the vortex core was initially about 10 percent of the 
rotor-blade chord, but more than doubled its size 
after one revolution of the rotor. According to vor- 
tex models that approximate the measured data, the 
core- radius circulation was about 79 percent of the 
large-radius circulation, and the large-radius circu- 
lation was about 67 percent of the maximum bound 
circulation on the rotor blade. On average, about 53 
percent of the maximum bound circulation resides 
within the vortex core during the first revolution of 
the rotor. 


* Presented at the American Helicopter Society 59th An- 
nual Forum, Phoenix, Arizona, May 6-8, 2003. Copyright 
©2003 by the American Helicopter Society International, Inc. 
All rights reserved. 
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ABBREVIATIONS 

CCD Charged coupled device 
CFD Computational fluid dynamics 
LDV Laser Doppler velocimetry 

PIV Particle image velocimetry 

INTRODUCTION 

A rotor wake follows a trajectory that repeat- 
edly passes near its own generating rotor blades be- 
fore finally being convected downstream. As a con- 
sequence, the rotor wake not only has a considerable 
influence on the performance of the rotor, but it is 
principally responsible for the high levels of noise and 
vibration that result from encounters between the ro- 
tor blades and their collective wake (ref. 1). Because 
of the complexity of the rotor wake, especially con- 
cerning the intensity and size of the dominant trail- 
ing vortex and its increasingly irregular trajectory, 
computational modeling of the flow field continues 
to be a challenge. Reliable experimental measure- 
ments are therefore needed to support such efforts. 

As new experimental techniques emerge, the ro- 
tor wake is frequently revisited in order to extract 
more detailed information about the structure of the 
trailing vortex (refs. 2 - 15). For example, not 
only does the stereo particle image velocimetry (PIV) 
measurement technique render all three components 
of velocity over a fairly large cross section of the flow, 
but the vector fields from a sequence of images taken 
during successive revolutions of the rotor can also be 
conditionally ensemble averaged to remove the ef- 
fects of vortex wander (ref. 14). Given this opportu- 
nity to acquire data more efficiently, highly resolved 
(in space and time) data on the development and ag- 
ing of the trailing vortex is now possible. In addition 
to enhancing our physical understanding of the rotor 
wake, specific characteristics (such as composition, 
size, intensity, and location) of the trailing vortex at 
different wake ages will enable computational fluid 
dynamics (CFD) models to be validated on a more 
detailed level. 

Many questions have been raised concerning un- 
expected features appearing within the trailing vor- 
tex (ref. 3), the composition of the vortex relative 
to blade loading (ref. 2) , vortex definitions (ref. 16), 
opinions dealing with the development and demise of 
vortex systems (ref. 17), and the appropriateness of 
simple vortex models (ref. 18). In response to these 
questions, an experiment was performed to provide 
data that focuses specifically on cross sections of the 
trailing vortex during a complete revolution of the 
rotor. The data are of sufficient quality to permit 
various modeling uncertainties to be resolved. This 


paper describes the test facility, model, flow mea- 
surement technique, and the data-acquisition and - 
processing procedure. A discussion of the results is 
provided along with specific findings. A considerably 
more detailed presentation is given in reference 19. 

TEST DESCRIPTION 

Model Rotor Setup 

A model rotor was placed in the center of the 
settling chamber of the 7- by 10-Foot Subsonic Wind 
Tunnel at Ames Research Center, and mounted hor- 
izontally with its axis of rotation coincident with the 
centerline of the test section. Aside from the influ- 
ence of the tunnel walls (approximately 4.6 m away) 
and the contraction into the test section (with an 
area ratio of 14:1), the rotor wake was allowed to 
pass through the test section without interference. 
A curtain was drawn across the air-exchanger sec- 
tion to prevent the wake from recirculating around 
the tunnel while fresh stagnant air was drawn by the 
rotor into the tunnel circuit (as required by mass 
continuity) . 

The rotor used in this experiment was a 1/7- 
scale model of the full-scale AH-l/OLS rotor (fig. 1). 
The rotor consisted of two blades, each with a chord 
of 103.9 mm, maximum thickness of 10.2 mm, and a 
span of 957.6 mm measured from the axis of rotation 
to the blade tip. Beginning at 25-percent radius, the 
twist of each blade was — 0.106°/centimeter (cm), so 
that the total variation in geometric twist along the 
span was 7.55° (from the center of rotation to the 
rotor tip). Collective pitch settings were referenced 
to the geometric pitch angle at 75-percent radius. 
Therefore, a collective pitch setting of zero meant 
that the blade pitch was about 5.0° at the root (held 
fixed for locations less than R/ 4) and —2.5° at the 
tip. The rotor hub had a built-in blade precone of 
1.30°. 

Stereo PIV Setup 

The primary elements of the PIV system con- 
sisted of a pulsed laser, optics to spread the beam 
into a thin sheet, charged coupled device (CCD) 
cameras to record the images, and a calibration tar- 
get. The plane of the light sheet was horizontal and 
passed just below the trailing edge of the rotor when 
the rotor was at zero azimuth (fig. 2). In order for 
both cameras to receive close to the same amount 
of forward-scattered light, two opposing, coplanar 
light sheets were established using synchronized 120 
millijoule (mJ) Nd:YAG lasers. Cameras and lasers 
were triggered off of a 1/rev reference signal that was 
keyed to the rotor. This signal was obtained by di- 
recting a small continuous-wave laser beam onto a 



small mirror attached to the hub of the rotor and 
intercepting the reflected beam with a fast-response 
photodetector. Acquisition software was used to de- 
lay the reference signal in order to obtain data at 
different wake ages. 

Images were acquired with 8-bit CCD cross- 
correlation cameras, each having a sensor array of 
2K x 2K pixels (nominal). Each pixel measured 7.4 
microns (/im) on a side. The cameras were operated 
in a double-exposure mode to acquire two noninter- 
laced, full-frame images during a single frame inter- 
val. The lenses were remotely translated to focus on 
the centerline of the image area, and the sensor (lo- 
cated inside the camera body) was rotated about its 
centerline to satisfy the Scheimpflug condition 1 . 

The lasers, sheet-forming optics, and cameras 
were rigidly mounted on a common structure, which 
could be traversed in three orthogonal directions. 
This arrangement allowed the cameras to be focused 
only on the trailing vortex region of the flow, and 
then traversed as required to follow the position of 
the vortex at different wake ages. 

Proper seeding of the flow is critical for obtain- 
ing accurate PIV measurements. The seed particles 
must be evenly distributed and of sufficient density 
to define the flow without altering its physical prop- 
erties. The particles must also be small enough to 
accurately follow the flow (especially challenging in 
accelerating flows), yet large enough so that they 
scatter a sufficient amount of light to be detected. 
The particle generator used in this test employed 
an inert gas to atomize a nontoxic, pharmaceutical- 
grade mineral oil. The mist was vaporized and then 
condensed before being released into the flow. The 
particle size was estimated by the manufacturer to 
be less than 0.5 //m. The particle generator was lo- 
cated upstream of the settling chamber so that the 
particles would mix with the air entering the tunnel 
circuit before being drawn into the wake of the rotor. 

PROCESSING WAKE IMAGES 

To obtain good resolution of the trailing vortex, 
a calculation grid containing 92 horizontal nodes by 
79 vertical nodes was constructed over the region of 
significant interest in the flow field. The area cov- 
ered by the grid resulted in approximately 12 pixels 
between nodes in both directions. The interroga- 
tion window was set at 22 pixels on a side, giving a 
45-percent overlap. The physical area measured in 
each image zone was 168 mm (with an axial incre- 
ment of 1.8 mm) by 106 mm (with a radial increment 
of 1.4 mm). The 1.8- by 1.4-mm increments, which 
defined the calculation step sizes in the axial and 

Uniform focus requires that the light sheet, camera lens, 
and camera sensor planes intersect along a common line. 


radial directions, also defined the resolution of the 
measurements. Therefore, the interrogation window 
(containing 484 pixels) covered an area measuring 
2.5 mm in the radial direction by 3.4 mm in the ax- 
ial direction. 

Consider an interrogation area in which 10 par- 
ticle pairs (the minimum number required to yield a 
valid resultant vector) have been registered on a 22- 
by 22-pixel image of a small region in a swirling flow. 
The software derives a representative in-plane veloc- 
ity vector based on a second-order calculation of all 
the individual particle displacements within that in- 
terrogation area. Note that the effective “probe size” 
is, therefore, less than the interrogation area. All the 
rotational velocity components within the vortex are 
accounted for, but they are represented in a collec- 
tive sense over contiguous regions defined by the size 
of the interrogation area. For an initial core diameter 
of 10 mm, the core would have been covered by about 
nine interrogation areas (78.5 mm 2 /8.5 mm 2 ). Con- 
sidering an interrogation overlap of 47 percent means 
that about 18 vectors would be calculated within the 
vortex core. Given that the size of the vortex grows 
as it ages, the number of vectors calculated within 
the core would also increase. In this present test, for 
example, a vortex that had grown to a diameter of 
20 mm would contain about 72 vectors. 

A single file for each wake age included a flag 
for every velocity measurement that indicated the 
nature of that vector (such as valid, invalid, interpo- 
lated, recalculated, or not calculated). Reduced data 
files containing the coordinates of each calculation 
node and the three components of velocity (in terms 
of displacement) were stored in ASCII format and 
then transferred to a mainframe computer for anal- 
ysis. The velocity components were first converted 
from displacement units to velocity units based on 
the pulse duration for that particular measurement. 
The physical coordinates were transformed so that 
r = z = 0 would correspond to the center of the 
rotor hub (fig. 2). 

Because of centripetal effects, heavier seed par- 
ticles appear to have been driven outward from the 
center of the vortex as it moved through the fluid. 
Because the size of the remaining lighter particles 
were most likely much smaller, less light would have 
been scattered from this region, thereby giving it the 
appearance of a “void.” Nevertheless, the particles 
that were present inside the void were sufficient to 
obtain PIV measurements. In searching for vortex 
characteristics that can be known with certainty, the 
size of the “apparent void” cannot always be equated 
to the size of the core, nor can the center of the 
void be assumed to be coincident with the center of 
swirl. On the other hand, the physical location of 


the void was an indisputable feature and was, there- 
fore, judged to be the most reliable indicator of how 
the vortex was convecting through the fluid. In addi- 
tion, in order to account for the ever-present vortex 
meander, the coordinates of the voids for a set of 
images acquired over time at a particular azimuth 
behind the rotor provided a very accurate basis for 
conditionally averaging the images. 

The distribution of particles around an apparent 
vortex void often exhibited a peculiar pattern. Parti- 
cles were not always uniformly distributed in the flow 
as it approached the rotor. In these instances, some 
form of particle “banding” frequently occurred adja- 
cent to the void region. Calculations have been made 
that show that banding is also a natural consequence 
of the balance of forces (drag, centrifugal, and Cori- 
olis) that act on the seed particles within a vortex 
(ref. 20) . Small particles settle in an annular region 
that is smaller than the vortex core, compared to 
larger particles that yield a band that is larger than 
the vortex core. Regardless of how uniform the flow 
appeared to be seeded in the present experiment, a 
fairly broad range of void sizes was observed at all 
wake ages (fig. 3). The average size of the void in- 
creased rapidly during the first 100° of azimuth, after 
which the rate of increase diminished substantially. 

In this study, the sequence for extracting the 
characteristics of the trailing vortex began with a 
physical examination of each raw image to deter- 
mine the size and location (relative to a superim- 
posed mesh) of the apparent void (fig. 4). At each 
wake age the images were averaged to smooth out 
the small irregularities that were invariably present 
in the flow. Because the vortex was located at dif- 
ferent positions from image to image (vortex mean- 
der), performing a simple average would have risked 
smearing out important details concerning the struc- 
ture of the vortex. Instead, the data were averaged 
only after aligning the set of images based on the 
most dominant recurring feature in the flow, which 
in this case was the apparent void. 

This procedure is referred to as conditional en- 
semble averaging, the condition here being the align- 
ment of all the images in the set based on their void 
centers. The image with a void center closest to the 
mean location that was established for the set of im- 
ages was selected as the “anchor.” The indices of 
all image matrices (that is, grid points i,j) for each 
wake age were then adjusted according to the offset 
of each void from the coordinates of the void in the 
anchor image. 

At this intermediate point, the conditionally av- 
eraged flow field at each wake age embodies the trail- 
ing vortex within a highly nonuniform downwash. 
The time-averaged downwash for a single-bladed ro- 


tor can be determined by averaging the flow field 
over an entire revolution of the rotor (refs. 5, 7). 
When the time-averaged downwash velocity field is 
subtracted from the phase- averaged velocity field, 
the resulting velocity profiles tend to appear more 
symmetric along the extremities of the vortex, but 
not near the core radius (ref. 5). In the present ex- 
periment, this procedure could not be implemented 
because data were not obtained over the required 
half revolution in each zone. However, the axial and 
radial locations of the voids at each wake age, along 
with the known time intervals between each wake 
age, were known and did yield the convection ve- 
locity of the vortex at each wake age. These com- 
ponents of the convection velocity, u c and w C} were 
subtracted from the phase- averaged flow field. How- 
ever, because the axial component of the downwash 
over the inboard region was significantly higher than 
the axial convection velocity of the vortex (u/ c ), ve- 
locity profiles across the center of the vortex cannot 
be expected to be symmetric. 

To extract information about the geometry of 
the vortex (after performing averages based on the 
void centers at each wake age), the center of the 
vortex was assumed to be at the center of swirl 2 . 
Using the void center as a starting point, the sur- 
rounding locations were interrogated to determine 
the best node for which the sum of the dot prod- 
ucts of two unit vectors produced a minimum over 
a neighborhood 3 of locations surrounding the candi- 
date node. The coordinates of the neighboring node 
relative to the candidate node defined one of the unit 
vectors. The velocity vector at the neighboring node 
defined the other unit vector. With the angles of 
these two unit vectors denoted by a and /?, the node 
nearest the center of swirl was said to be found when 
the following expression was a minimum: 

N 

£|cos(a-/?)|/tf (1) 

n—l 

where N is the total number of neighboring points 
considered. Relative to the center of swirl for a pure 
vortex in a stationary flow, a and /? will be orthog- 
onal for all neighboring points and the sum will be 
exactly zero. For most real flows the sum will not be 
precisely zero, and the procedure works best when 
the neighborhood under consideration does not ex- 
tend beyond about one or two core diameters (which 
is normally sufficient to cover most of the vorticity 

2 Neither the void center nor the swirl center was necessarily 
at the center of the vorticity that was accreted from the rotor 
wake, itself having the form of a spiraling sheet containing 
varying amounts of vorticity. 

3 Because the quality of the data within the void may be in 
question, these vectors were omitted during this operation. 



in a single trailing vortex). 

The size of the vortex core was found by first di- 
viding the surrounding region into annular zones and 
calculating the average swirl velocity and the associ- 
ated mean radius of each annulus. These data were 
fit with a least-squares spline subject to the condi- 
tion that the resulting curve have only one inflection 
along its inner extent (nominally set to twice the esti- 
mated core radius). The point of inflection was used 
to define the radius of the vortex core, because that 
is where a maximum value for the swirl velocity is 
reached. This procedure was considered to offer a 
more rational approach for determining the size of 
the vortex, especially when the velocity peaks that 
are characteristic of vortices (and upon which core 
sizes are traditionally based) are dependent on how 
the vortex is sliced (vertical, horizontal, or other- 
wise) . 

Vorticity was calculated by integration (circula- 
tion box method) according to: 

u >0 ~ A~ 1 T = A' 1 <j)(u } w) ■ d(f } z) (2) 

where the direction of integration was such that the 
enclosed area, A , was on the left of the integration 
path. The size of each circulation box was defined 
by the local grid spacing, which on average rendered 
A = 2.5 mm 2 . 

TEST CONDITIONS 

The experiment was performed at a constant ro- 
tor speed of 870 rpm (14.5 Hz), which for a rotor di- 
ameter of 191.5 cm corresponds to a tip speed of 87.2 
m/s. Based on an average ambient temperature of 
16.7 degrees Centigrade (°C), the Reynolds number 
based on rotor-tip chord was R e = 0.6 x 10 6 and the 
Mach number was M = 0.26. The collective pitch 
angle of the rotor was fixed at 11°. 

To obtain high-resolution PIV images, the cam- 
eras were focused on a small region of the wake be- 
hind the trailing edge of the rotor. The camera view 
was radially and axially translated to provide three 
zones of coverage (A, E, and F). The laser pulse de- 
lays needed to capture the highest velocities ranged 
from 3.5 to 5.0 /is. Results presented later in this re- 
port show that even for the longest pulse delay, 5 mi- 
croseconds (fj, s), the vortex moved less than 1/1000 of 
the core radius between pulses. The blade was con- 
sidered to be at zero azimuth when the trailing edge 
of the rotor blade was horizontal. When the blade is 
in this position, the light sheets (both of which in- 
tersect in a horizontal plane) pass just beneath the 
trailing edge of the blade. The procedure for taking 
images at the different wake ages was accomplished 


by delaying the 1/rev signal from the rotor, thereby 
allowing the tip of the rotor blade to move past the 
measurement plane by the desired azimuthal incre- 
ment. The data were taken at a nominal increment 
of 15° (fig. 5), with some exceptions as required to 
avoid reflections of the laser light from any portion 
of the blade surface that might damage the camera 
sensors. 

Rotor loads and tunnel velocity measurements 
(table 1) were taken from an earlier experiment (ref. 
15). Strictly speaking, the rotor is not in hover, as 
evidenced by the magnitude of the velocity upstream 
of the rotor. Because the flow around the rotor is 
bounded by the wind tunnel walls, and because the 
rotor wake lowers the pressure in the test section (in 
turn inducing a finite flow everywhere ahead of the 
rotor), a slight climb condition cannot be avoided 
(see ref. 15 for a full discussion). Assuming the up- 
stream velocity entirely represents the climb velocity, 
then the climb rate for this part of the experiment 
can be calculated from Woo/RQ. 

Table 1: Rotor loads 


Rotational speed of rotor, rpm 

870 

Collective pitch angle, deg 

11 

Solidity 

0.0691 

Torque, foot-pounds (ft-lbs) 

11.69 

Thrust, lb 

39.76 

Upstream velocity, ft/s 

1.28 

Climb rate 

0.0045 

Hover figure of merit 

0.61 

Ct 

0.00657 

Cq 

0.00062 

Ct/g 

0.095 


RESULTS AND DISCUSSION 

The center of the vortex void, the center of vor- 
ticity, the center of swirl, and the center between 
velocity peaks are all considered to be distinct quan- 
tities. For each wake age, starting with = 0°, the 
flow field has been summarized in terms of quantities 
that focus on the trailing vortex (i.e., fig. 6). Vortex 
wander was qualitatively identified by plotting the 
locations of the vortex voids for each image set (up- 
per left comer of figure). To the right of this plot is a 
sketch of the rotor showing the location of the rectan- 
gular zone covered by the entire calculation grid (92 
steps in the z direction and 79 steps in the r direc- 
tion). The three field plots in the left-most column 
show the out-of-plane vorticity, Qq (contours with red 
and blue hues to indicate direction and a continuous 
range of saturations to indicate magnitude), inplane 



velocity relative to the convection velocity of the vor- 
tex, u r + (vectors with a continuous range of hues 
to indicate magnitude), and the out-of-plane veloc- 
ity, v (contours with red and blue hues to indicate 
direction and a continuous range of saturations to 
indicate magnitude). The three profiles in the cen- 
ter column show w r , and v along z = Constant 
and passing through the center of swirl. The central 
portion of these profiles is shown in red and demar- 
cated with dashed vertical lines in order to identify 
the void region. The horizontal and vertical dashed 
lines included in each of the field and profile plots 
mark the location of the swirl center. 

The lower plot in the right-most column shows 
the inplane swirl velocity, <p t averaged over incremen- 
tal annular regions with radii (C) measured from the 
swirl center. The upper plot in the right-most col- 
umn shows an enlargement of the inplane velocity 
(u r -f w r ) } limited to a region twice the diameter of 
the vortex core. The two circles identify the size 
of the vortex core (solid curve) and the vortex void 
(dashed curve). Finally, values are given at each 
wake age for the void center, the vorticity center 
(or centroid of vorticity), the swirl center, the vor- 
tex convection velocity (obtained from the path of 
the vortex voids), and the core radius along with the 
maximum swirl velocity (both of which are fixed by 
the peak value found from the curve fit to <p vs. £)• 
These values are also given in tables 2 to 4. 

Based on a review of the data for all wake ages 
(selected cases are given in figs. 6-10), the highest 
velocities are found to occur within the first 15 de- 
grees behind the rotor blade. During this interval, 
the maximum inplane (swirl) and out-of-plane ve- 
locities reached 44 and 12 percent of the rotor tip 
velocity (87.2 m/s), respectively. Although the de- 
gree of vortex wander was somewhat erratic, it gen- 
erally increased with wake age. On average, the 
standard deviation (based on the mean location of 
the void center) was only about 9 percent of the di- 
ameter of the vortex core (the worst case being 18 
percent). The calculated locations of the vorticity, 
swirl, and velocity-peak 4 centers generally differed 
from the void center by about 11, 26, and 11 percent 
of the vortex diameter, respectively (fig. 12). There 
was no apparent correlation between the void and 
swirl centers from image to image. 

The out-of-plane velocity exhibited the most un- 
usual appearance. Aside from the boundary layer 
wake that is shed from the trailing edge of the ro- 
tor blade and is most distinct when xp < 30° , there 
appear to be two concentrated zones of oppositely 
directed flow near the center of the vortex. Because 

4 The velocity-peak center is based on the locations that 
define d h (see fig. 11). 


both of these zones are generally offset from the vor- 
tex center, the bisecting profiles taken along z = 
Constant do not fully reveal the magnitude of the 
“excess” and “deficit” velocities that contribute to 
the overall structure of the trailing vortex. 

Even though the convection velocity of the vor- 
tex at each wake age has been subtracted, the rela- 
tive “vertical” component of velocity, w r , exhibits an 
asymmetric profile. Greater symmetry near the vor- 
tex may be achieved by subtracting the local down- 
wash, which can only be approximated from a time 
average of the flow field at corresponding locations 
over one or more revolutions (refs. 5, 7). However, 
it is not clear that this effort will yield any better 
understanding of the vortex structure, especially be- 
cause the resulting profiles (including the peaks) will 
still be asymmetric (refs. 5, 10). Given the inevitable 
likelihood of asymmetry, traditional estimates of the 
core size and peak velocity would still have to be 
based on an average between the two halves of the 
asymmetric profile, and this is effectively what is al- 
ready provided in the more general annular averages 
that produce the <p vs. C plot. Nevertheless, the pro- 
file plots are important because they retain certain 
nuances that can prove helpful in explaining other 
unexpected features in the flow. 

During the initial development of the trailing 
vortex, ip < 30°, the vorticity within the vortex 
reached a peak value very near the center of the vor- 
tex, and rapidly diminished to near zero at a distance 
not far beyond the core diameter. For larger wake 
ages, V* > 45°, the peak values of vorticity appear to 
reside in an annular region defined by the size of the 
void, yet are clearly within the core diameter. As- 
suming these measurements are valid, this variation 
in vorticity would imply that much of the interior 
region of the vortex was not fully in solid-body rota- 
tion. Evidence for this can be found in the vertical 
component of velocity, w ry which shows a change in 
the slope near the vortex center (see, for example, 
the profile at ip = 195° appearing in fig. 9). This 
peculiar departure from the traditional view of a vor- 
tex has appeared in other rotor vortex data (ref. 6, 
where the flow was seeded with incense particles and 
measurements made with a laser Doppler velocime- 
try (LDV) technique, and ref. 24). 

To be sure, other explanations (seeding being 
among them) for the unusual velocity and vorticity 
profiles over the interior of the vortex are plausible. 
Although centripetal effects may have been a factor, 
the extent to which the velocity within the core was 
significantly altered is unknown because the size of 
the particles was not measured. Another factor may 
have been the rival between legitimate, but barely 
visible, small particles and the presence of sensor 



noise, the result being an abundance of interpolated 
vectors within the void region. Until these issues 
can be more thoroughly addressed, it is uncertain 
whether the data show an onset of vortex instability; 
reflect a difference in particle dynamics because of a 
nonuniform distribution of particle sizes in a highly 
accelerated flow; or are simply due to a processing 
limitation arising from a marginal scattering of light 
from very small particles that resulted in incorrect 
interpolations. 

The location with the greatest number of in- 
terpolations always occurred within the void region. 
The number of interpolated values contributing to 
the average at these “worst-case” locations ranged 
from 16 to 60 percent. Although far from ideal, no 
location consisted solely of interpolated values. 

To address a possible flaw in the logic used to 
average the data during the analysis phase, the case 
for ip = 195° was reexamined. First, the results for 
single images (hence, unaveraged) were individually 
studied to see if the nonlinear slope in the core ve- 
locity might simply have been an artifact of the av- 
eraging process. Although the un averaged velocity 
profiles were noticeably more “noisy,” the peculiar 
change in slope near the center of the vortex was 
found to be present and similar to the conditionally 
averaged results. As a final check, a different cri- 
terion was considered for obtaining the conditional 
average. In this case the more classical approach of 
using the midpoint between the velocity peaks (high- 
est maximum and minimum values occurring dur- 
ing horizontal cuts across the vortex) was adopted. 
Again, the variation in slope across the vortex was 
essentially unchanged from that observed when the 
conditional averaging criterion was based on the lo- 
cations of the voids. In fact, the same conclusion was 
reached after considering the swirl centers as the av- 
eraging criterion. 

Another curious feature of the vortex is its 
slightly elliptical shape. This was also observed in 
the flow pattern constructed from LDV measure- 
ments in the near wake of a hovering rotor (ref. 22) A 
minor amount of ellipticity can be attributed to the 
fact that the axis of the vortex is not perpendicular 
to the plane of the light sheet. This would cause the 
vortex to appear stretched in the radial and axial 
directions (relative to the rotor coordinates). The 
amount of stretching in both directions was nearly 
the same during the early wake ages, thereby in- 
creasing the apparent core size by about 0.06 per- 
cent and having essentially no effect on ellipticity. 
As the vortex aged, the apparent stretching became 
progressively more biased in the axial direction (be- 
cause the radial displacement was approaching a con- 
stant), but the degree of ellipticity caused by the an- 


gle of the vortex axis was still less than 0.4 percent 
after one revolution. Judging by the amount of ellip- 
ticity observed in the data (sometimes reaching an 
aspect ratio of about 1.23), and by the alternating 
orientation of the major axis, the cause for this ellip- 
tical flow pattern cannot be due to a changing angle 
of the vortex axis as the vortex ages (unless small- 
amplitude undulations are superimposed on the vor- 
tex trajectory, although the data in figure 12 do not 
show this). 

To study this phenomenon in more detail, the 
dimensions of the vortex (major and minor axes) 
were equated to the distances between the velocity 
peaks found from horizontal and vertical cuts across 
the center of the vortex (recall fig. 11). The peak- 
to-peak distances associated with the horizontal and 
vertical cuts were designated by dh and d Vy respec- 
tively, and these are shown for wake ages ip — 0° to 
390°. Using red and green color bars to indicate the 
size of dh relative to d V) the data suggest that the 
elliptical pattern is rotating. This figure also shows 
that the nonlinear behavior of the velocity profile 
within the core does not depend on the direction of 
the cut (horizontal or vertical) or the direction of the 
major axis of the ellipse (horizontal or vertical). 

According to the location of the void center, the 
radial and axial coordinates of the vortex with in- 
creasing wake age are known (fig. 13). Based on 
curve fits to this data, the corresponding convection 
velocities were calculated. These data indicate that 
the highest radial velocity occurred immediately af- 
ter the trailing vortex left the trailing edge of the ro- 
tor. This figure shows that the magnitude of the ra- 
dial velocity rapidly decreased and then approached 
an asymptotically diminishing value near the com- 
pletion of one rotor revolution. The behavior of the 
axial velocity component was quite different. The 
magnitude of the axial velocity was lowest imme- 
diately after the trailing vortex left the rotor, but 
then increased until its progress was noticeably im- 
peded after the passage of the following rotor blade. 
The velocities may have a slightly different appear- 
ance if curve fits to the data are separated into seg- 
ments before and after passage of the second blade 
(rp = 180°). The wake contraction can also be viewed 
in the r-z plane (fig. 12), which shows the character- 
istic “knee” due to the influence of the second blade, 
and the increasing dominance of the axial velocity as 
the wake aged. 

Based on the radial distribution of the swirl ve- 
locity at each wake age, the peak (or inflection) in 
the curve fit was used to define the core size and 
the maximum swirl velocity. The initial diameter of 
the vortex core was about 10 percent of the chord of 
the rotor blade (c — 104 mm) and more than dou- 



bled its size after one revolution of the rotor (fig. 
12). The maximum swirl velocity was highest during 
the initial formation of the trailing vortex, decreased 
rapidly during the first half revolution, and settled 
on a nearly constant value (on average) during the 
last half revolution of the rotor. It is difficult to say 
exactly how the intensity of the vortex (which is re- 
lated to the maximum swirl velocity) changed with 
wake age because of the unavoidable influence of the 
downwash on the magnitude of the swirl velocity (the 
vortex is embedded in a nonuniform velocity field). 

To examine more closely the structure of the 
trailing vortex, contour plots of the vertical compo- 
nent of velocity, w r> and the out-of-plane vorticity, 
J 0) are shown at wake ages ip = 15° and 195° along 
with profiles taken from several closely spaced hori- 
zontal cuts across the vortex (figs. 14). The velocity 
profiles are clearly more linear within the core of the 
vortex at ip = 15° than they are at ip = 195°. The 
corresponding vorticity profiles exhibit a single peak 
when crossing the vortex at ip = 15°, whereas a dou- 
ble peak is evident along several cross sections of the 
flow at ip = 195°. 

Another point of interest is the radial distribu- 
tion of circulation across the vortex. Using the swirl 
velocity to calculate the circulation, and normalizing 
both (p and T by values existing at the core radius 
(Co), results for ip — 15°, 195°, and 375° are shown in 
figure 15. As the distance increased from the center 
of swirl (C = 0), the circulation rapidly increased un- 
til reaching a point just beyond the core radius, after 
which the rate of growth decreased significantly. 

Similar to the selected views of circulation, the 
radial distribution of vorticity was examined during 
the development of the trailing vortex. However, be- 
cause the vorticity center and the swirl center (as 
well as the void center) data were generally not collo- 
cated, the normalized vorticity profile was referenced 
to the vorticity center. The peak vorticity occurred 
at the center ( k = 0) when ip = 15° (fig. 16). At 
subsequent wake ages (ip = 195° and 375°), the vor- 
ticity in the center decreased and the peaks progres- 
sively shifted toward the core diameter. A similar 
change in the distribution of vorticity was reported 
in an earlier fixed-wing experimental study using a 
vorticity meter (ref. 23). As discussed earlier, it is 
not known to what extent this characteristic might 
be due to potentially inaccurate data within the void 
region. Nevertheless, the fact that this radial shift 
in the vorticity peak has been observed while using 
other measurement techniques (vorticity meter and 
LDV) suggests that the PIV data are authentic and 
that the vortex was indeed becoming unstable. 

The self-similar nature of the swirl velocity and 
the vorticity for the entire set of wake ages (ip = 0° 


to 390°) are shown in figure 17. The swirl veloc- 
ity was made dimensionless by its value at the core 
radius ( <p/<p c ), and the radius (measured from the 
swirl center) was made dimensionless by the core ra- 
dius (C/Cc)- Some data points for the dimensionless 
swirl, Tp, exceed 1.0 because the maximum swirl ve- 
locity, p>c, was based on curve fits to <p vs. £ at each 
wake age. The dimensionless form of a popular vor- 
tex model (ref. 18) 



has also been included in the plot. The best fit to the 
data occurred when n = 3 and r*/T c = 1.26. Ac- 
cording to this model, the circulation around a con- 
tour placed at the core boundary, r c , was about 79 
percent of the total far-field (or large-radius) value, 
r*. Note that at the core radius Tp = £ = 1, so that 
equation (3) requires that F*/r c = 2 L / n = 1.26 when 
n = 3. 

The self-similar distribution of accumulated vor- 
ticity can be found from the product of the local 
vorticity (w$) and the incremental area over which 
it acts (ArAz), and then summed over regions de- 
fined by the radius from the vorticity center («). In 
principle, the result should be equivalent to the vari- 
ation of circulation with radius, provided the swirl 
and vorticity centers are coincident. Within the core 
( k < Cc), the distribution of vorticity is quite similar 
and is well represented by the model, TpQ. Beyond 
the core radius the accumulated vorticity asymp- 
totically approaches a constant value and on aver- 
age agrees with that predicted by the vortex model, 
which states that — >• 1.26 as C becomes large. 

The core-radius circulation, r c , was observed to 
increase slightly with wake age (fig. 18). Because 
some vorticity was found to exist beyond the core 
radius at all wake ages, and because the core ra- 
dius generally increased with wake age (fig. 18), it 
is probable that additional vorticity was gradually 
being gathered into the core. The total far-field vor- 
ticity appears to be fairly constant, T* = 1.42 m 2 /s. 
The growth in the size of the vortex core can be 
approximated by the following model (ref. 24, and 
further developed in ref. 22): 

Cc = (Coo + 5(^/0) 1/2 (4) 

where 8 = 16 provides the best fit to the present 
data (for laminar flow, 8 = 1). Data included from 
another experiment (ref. 22, based on LDV measure- 
ments) was reported to be best fit by the model when 
8 = 8. It should be noted that <5 strongly depends 
on the initial core size, CcO, and that lower values of 



Cco require lower values of S. The results presented 
in reference 25, which are for a single-bladed rotor, 
show a similar range of core sizes and growth rate. 

Based on an isolated vortex model (refs. 20, 
25) that assumes a turbulence mixing length that 
is proportional to the vortex radius, the ratio of the 
core-radius to large-radius circulation should be fixed 
by the vortex Reynolds number, T+/i/. However, the 
data for this test {T^/u » 1.2 x 10 5 ) fall above the 
curve and are, on average, aligned with a laminar 
value for r c /r* (fig. 18). This would imply that 
S & 1 rather than 16 as required for a good fit to 
the Cc V' data, and would suggest that either the 
models are imprecise or that the core radius is not 
as large as measured. 

Just as the core-radius circulation was less than 
the large-radius circulation during the first revolu- 
tion of the rotor, the large-radius circulation is ex- 
pected to be less than the maximum bound circu- 
lation on the blade. Although the blade circulation 
was not measured in this test, a value can be ap- 
proximated from blade-element momentum theory 
(ref. 26). For a blade with constant chord, the to- 
tal thrust coefficient can be stated in terms of the 
dimensionless inflow, A, as 

C T = f 4A 2 rdr (5) 

Jo 


expressions for the inflow in the form of A = A 0 r n 
have been proposed (refs. 20, 26). The results for 
n — 0 (constant inflow), n = 0.5 (linear circulation), 
and n = 1 (linear inflow) are shown in figure 18. 
Although this power function has served well to il- 
lustrate various performance trends, the inflow used 
in the present analysis was 


A = A 0 ysin~r (9) 

From equations (5), (7), and (9) the maximum bound 
circulation on the rotor blade is 


(IL) 

\ClRe) t 



(10) 


Recalling that Ct = 0.00657 and a = 0.0691, then 
equation (10) states that (r^ma* = 2.13 m 2 /s for 
this test. To summarize: 


~ = 0-79 (11) 

= 0-67 (12) 

Lbmaz : 

implying that the core-radius circulation is only 
about 53 percent of the maximum bound circulation 
on the rotor blade. 


and the sectional value in terms of the local circula- 
tion, r 6 , on the blade as 

JC - = m riT (6 > 

From the above two equations it follows that the 
inflow velocity and the blade circulation (both are 
functions of r) are related by 


o-T b 

QRc 


= 4A 2 


( 7 ) 


The local inflow velocity for hover can be stated as 
(ref. 26) 



where 9 is the local blade angle. Because the climb 
rate in this test was insignificant, hover conditions 
are assumed to apply. The change in blade angle 
along the span for the rotor used in this test is shown 
in figure 18. If the blade angle, <9(r), is known, and if 
it is assumed that Ci a ^ 0.1, the local inflow veloc- 
ity predicted by equation (8) can be determined (fig. 
18). However, in order to more easily perform the 
integration given in equation (5), simple analytical 


CONCLUSIONS 

1 . The highest velocities occurred within the first 
15° behind the rotor blade. The swirl and axial ve- 
locities reached 44 and 12 percent of the rotor tip 
speed, respectively. 

2. The initial diameter of the vortex core was 
about 10 percent of the blade chord and more than 
doubled its size after one revolution of the rotor. 
The rate of growth followed a square-root function 
of wake age. 

3. The boundary , layer wake that was shed along 
the span of the rotor blade disappeared after t/> = 
30°; however, two concentrated zones of oppositely 
directed axial flow jvere observed near the center of 
the vortex throughout the first rotor revolution. 

4. The void, vorticity (or centroid of vorticity), 
and the swirl centers occurred at distinct locations, 
but were never separated by more than 26 percent of 
the vortex diameter. 

5. The axially directed velocity profile, relative to 
the convection velocity of the void center, was asym- 
metric because of the influence of the rotor down- 
wash. The vortex also appeared to have an elliptical 
shape that rotated as the wake aged. 


6. The peak vorticity occurred near the center of 
the vortex during the early development of the trail- 
ing vortex, and then progressively shifted toward the 
core diameter. It is presently unknown to what ex- 
tent this result is due to weakly defined characteris- 
tics within the void region. 

7. The magnitude of the radial distribution of the 
circulation rapidly increased until reaching a point 
just beyond the core radius, after which the rate of 
growth decreased significantly. The core-radius cir- 
culation increased slightly with wake age, but the 
large-radius circulation appeared to remain relatively 
constant. 

8. The radial distributions of swirl velocity and 
vorticity exhibit self-similar behaviors, especially 
within the core. According to vortex models that 
approximate the measured data, the core-radius cir- 
culation was about 79 percent of the large-radius cir- 
culation, and the large-radius circulation was about 
67 percent of the maximum bound circulation that 
was calculated for the rotor blade. On average, about 
53 percent of the maximum bound circulation resides 
within the vortex core during the first revolution of 
the rotor. 
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Table 2: Locations related to vortex 


Wake age 
(deg) 

Vorticity center 
r (mm) z (mm) 

Swirl center 
r (mm) z (mm) 

Void center 
r (mm) z (mm) 

Core radius 
(mm) 

Void radius 
(mm) 

0 

945.78 

2.89 

946.11 

1.93 

946.00 

3.00 

5.00 

0.71 

15 

934.78 

-1.01 

934.25 

-1.31 

934.70 

-0.80 

5.25 

2.21 

30 

922.28 

-7.88 

922.99 

-8.66 

922.10 

-8.80 

5.00 

2.85 

45 

909.53 

-16.08 

909.51 

-14.88 

910.00 

-16.00 

6.75 

3.42 

60 

899.23 

-26.00 

899.35 

-25.80 

899.30 

-26.70 

6.75 

3.74 

75 

888.40 

-35.37 

888.98 

-34.31 

888.40 

-35.20 

6.75 

4.42 

90 

877.85 

-44.60 

878.04 

-45.51 

878.20 

-45.20 

7.50 

5.04 

105 

868.99 

-54.56 

867.06 

-54.17 

869.00 

-55.20 

7.75 

5.44 

120 

858.57 

-64.35 

858.66 

-63.02 

858.50 

-64.50 

8.50 

5.55 

135 

850.34 

-76.28 

851.04 

-73.69 

850.00 

-75.70 

8.50 

5.84 

150 

842.34 

-83.49 

839.95 

-84.16 

842.20 

-85.40 

10.50 

5.94 

180 

829.99 

-112.12 

828.45 

-112.83 

830.60 

-111.80 

8.75 

6.71 

195 

820.94 

-130.12 

820.29 

-130.54 

823.60 

-131.10 

8.75 

7.08 

210 

820.87 

-148.14 

819.50 

-146.08 

821.10 

-149.40 

9.51 

6.80 

225 

814.53 

-165.86 

812.56 

-166.17 

816.60 

-166.00 

8.75 

6.82 

240 

810.73 

-185.65 

808.70 

-183.76 

812.70 

-185.30 

9.99 

6.62 

255 

809.76 

-200.91 

808.47 

-201.46 

810.50 

-202.00 

10.24 

7.49 

270 

810.10 

-218.07 

807.03 

-217.92 

810.70 

-217.60 

9.75 

6.51 

285 

804.78 

-231.66 

801.58 

-231.76 

804.50 

-232.60 

10.00 

6.28 

300 

802.36 

-246.28 

800.52 

-245.89 

803.10 

-247.00 

9.75 

6.57 

315 

802.14 

-262.10 

800.28 

-261.70 

802.40 

-263.10 

9.75 

6.72 

328 

801.18 

-277.29 

803.30 

-276.32 

801.10 

-277.00 

10.00 

7.04 

344 

794.13 

-293.91 

791.73 

-295.04 

794.50 

-293.60 

10.50 

6.76 

375 

795.78 

-312.19 

796.70 

-309.02 

795.50 

-311.90 

11.00 

7.11 

390 

798.17 

-328.52 

797.79 

-325.25 

798.90 

-329.10 

11.25 

7.14 



Table 3: Vortex convection and shape 


Wake age 

(de«) 

Vortex convection 
u c (m/s) xjj c (m/s) 

Max swirl velocity 
(m/s) 

Velocity peak separation 
d h (mm) d v (mm) 

0 

-3.43 

-3.03 

30.77 

10.77 

9.07 

15 

-3.43 

-3.04 

35.36 

9.58 

11.15 

30 

-3.41 

-3.07 

32.24 

10.93 

11.02 

45 

-3.37 

-3.14 

24.88 

13.64 

12.88 

60 

-3.24 

-3.31 

25.86 

13.51 

14.55 

75 

-3.23 

-3.39 

27.39 

13.68 

12.75 

90 

-3.12 

-3.57 

24.59 

13.65 

14.58 

105 

-2.99 

-3.79 

23.26 

13.67 

14.56 

120 

-2.84 

-4.04 

20.29 

13.49 

18.30 

135 

-2.66 

-4.30 

16.62 

14.91 

18.36 

150 

-2.48 

-4.57 

17.07 

16.26 

18.19 

180 

-2.09 

-5.07 

20.78 

17.61 

14.46 

195 

-1.89 

-5.28 

17.50 

19.03 

18.14 

210 

-1.70 

-5.43 

17.94 

18.04 

17.81 

225 

-1.52 

-5.53 

17.52 

17.99 

15.68 

240 

-1.35 

-5.58 

15.32 

19.37 

17.65 

255 

-1.20 

-5.58 

17.82 

19.38 

19.71 

270 

-1.06 

-5.53 

17.27 

19.02 

18.26 

285 

-0.95 

-5.45 

15.94 

18.99 

16.44 

300 

-0.85 

-5.35 

19.57 

18.97 

18.32 

315 

-0.77 

-5.25 

18.20 

19.04 

18.10 

328 

-0.72 

-5.17 

17.71 

19.07 

18.26 

344 

-0.68 

-5.08 

17.84 

19.02 

18.29 

375 

-0.64 

-4.99 

16.35 

18.96 

21.74 

390 

-0.64 

-4.99 

15.63 

19.11 

21.82 


Table 4: Vorticity and velocity extrema 


Vorticifcy normal to Velocity within Velocity normal to 


Wake age measurement plane measurement plane measurement plane 


(deg) 

max (1/s) 

min (1/s) 

max (m/s) 

min (m/s) 

max (m/s) 

min (m/s] 

0 

20463.3 

-2484.1 

36.32 

1.08 

4.79 

-10.78 

15 

22942.8 

-994.0 

38.76 

1.21 

2.62 

-9.40 

30 

21014.8 

-1290.1 

35.08 

1.54 

1.83 

-9.25 

45 

11374.9 

-1198.1 

26.59 

0.23 

2.33 

-7.61 

60 

11808.0 

-1055.6 

26.67 

0.02 

2.79 

-7.16 

75 

11929.1 

-1043.4 

29.03 

0.14 

2.71 

-6.72 

90 

10505.9 

-918.7 

26.52 

0.05 

4.57 

-6.61 

105 

9416.3 

-1028.0 

25.78 

0.05 

4.49 

-6.56 

120 

7657.5 

-1050.3 

23.93 

0.03 

4.69 

-6.55 

135 

5812.4 

-1172.7 

22.56 

0.20 

5.03 

-6.96 

150 

6990.6 

-1288.0 

25.69 

0.85 

4.55 

-9.30 

180 

8575.8 

-1143.9 

24.23 

1.37 

3.31 

-5.31 

195 

8755.1 

-1814.7 

23.55 

0.72 

3.52 

-6.09 

210 

8154.6 

-900.5 

23.14 

0.61 

3.48 

-5.76 

225 

7746.3 

-1147.3 

22.40 

0.44 

3.26 

-4.50 

240 

6852.9 

-1623.4 

19.61 

0.24 

3.46 

-4.41 

255 

7020,4 

-724.5 

20.62 

0.76 

4.19 

-5.36 

270 

7006.7 

-870.7 

22.13 

0.27 

2.77 

-4.82 

285 

6732.1 

-1318.6 

19.52 

0.06 

2.40 

-4.43 

300 

7062.0 

-706.5 

22.23 

1.37 

2.11 

-5.33 

315 

6382.8 

-773.3 

19,93 

1.20 

2.35 

-3.68 

328 

6532.7 

-788.4 

21.20 

0.90 

2.79 

-4.44 

344 

6841.8 

-859.0 

20.49 

0.38 

1.72 

-9.15 

375 

5851.2 

-492.8 

.20.95 

0.73 

2.63 

-7.30 

390 

5977.2 

-734.8 

21.39 

0.05 

4.32 

-7.26 
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Figure 1: OLS model rotor blade. 



Figure 2: Setup for PIV measurements. 
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Figure 5: Wake ages at which PIV data were recorded in zones A, E, and F. 
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Figure 6: Velocity and vorticity characteristics with reference to the center of swirl at jp = 0°. 
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Figure 7: Velocity and vorticity characteristics with reference to the center of swirl at ^ = 30°. 
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Figure 10: Velocity and vorticity characteristics with reference to the center of swirl at tft = 375°. 



Figure 11: Evidence for the dynamic elliptic shape of the vortex during ip = 0° to 390 
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Figure 12: History of vortex attributes from tb = 0° to 390°. 







Figure 13: Axial and radial convection of the trailing vortex from ^ = 0° to 390 
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Figure 14: Velocity and vorticity profiles along several cross sections of the trailing vortex at V = 15° and 
195°. 
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Figure 16: Vorticity (w«), calculated along radius /c relative to the vorticity center using the inplane velocity 
field at V- = 15°, 195°, and 375°. 
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Figure 18: Vortex size and strength during wake ages ip = 0° to 390° and hover inflow models. 



